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We review some surprising links which have been discovered in the last few years 
between the theory of certain ordinary differential equations, and particular inte- 
grable lattice models and quantum field theories in two dimensions. An application 
of this correspondence to a problem in non-Hermitian (PT-symmetric) quantum 
mechanics is also discussed. 
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1. Introduction 

In the last few years, detailed connections have begun to emerge between two previously- 
separated areas of mathematical physics: the spectral properties of ordinary differential 
equations such as the Schrodinger equation, and the study of integrable lattice models and 
integrable quantum field theories using techniques related to the Bethe ansatz. This is 
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sometimes given the (perhaps over-grand) title of the 'ODE/IM correspondence'. While the 

correspondence remains at the level of a mathematical coincidence, albeit a detailed one, it 
has already provided some important new insights into a variety of problems. 

On the 'ODE' (ordinary differential equations) side, some of the key themes are Her- 
mitian and non-Hermitian spectral problems, Schrodinger equations, the WKB method, 
analytic continuation and T'T-symmetry. These have been developed by a variety of au- 
thors; some important names and references for the story that we want to tell are Sibuya [1] , 
Voros [2], Bessis and Zinn- Justin [3], and Bender and Boettcher [4,5]. 

The 'IM' (integrable models) aspect has a similarly long history. The theories relevant 
for the first example of the correspondence to be discovered [6] are the six-vertex model, 
and the quantum field theory of the massless sine-Gordon model. The Bethe ansatz, TQ 
relations and fusion hierarchies all have a role to play, and the works on these topics by 
Baxter [7], Kliimper and Pearce [8], and Bazhanov, Lukyanov and Zamolodchikov [9-11] 
are particularly relevant. 

In a short survey such as this, it is impossible to do justice to all of this background. In 
the nc;xt two sections some of the basic vocabulary will be sketchc;d. More dc;tailed reviews 
have been given in [12,13], and an extended version is currently in preparation. Subsequent 
developments following the initial observation in [6] can be found in [14-28]. 

2. First prologue: the basics of integrable models 

Consider a rectangular two-dimensional lattice, of size N x M. On each link of the lattice, 

attach an arrow, or spin, a, pointing in one of the two possible directions along the link. 
An assignment of a spin to each link i of the lattice gives a configuration {ui}. If we 
impose the additional constraint that the number of arrows pointing into each vertex is 
equal to the number pointing out, then we are on the way to defining the six-vertex model. 
The main task is to calculate averages over all possible configurations, weighted by relative 
probabilities which can be found as follows. First, the local Boltzmann weights W[a,b,c,d] 
must be specified, one number for each possible set of values taken by the spins on the 
four links connected to a single site of the lattice. Then the relative probability of a given 
configuration is simply the product of these local weights over all sites. A key quantity is 
now the sum of these relative probabilities over all possible configurations, which gives the 
normalisation factor for the calculation of the actual probability of any quantity. Otherwise 
known as the partition function, it is given schematically as 



The local Boltzmann weights W, and hence also the partition function Z, may depend on 
physical parameters such as the temperature. For the six-vertex model there is just one 
physically-relevant quantity to specify, namely the anisotropy rj. In the related quantum 
field theories this becomes the sine-Gordon coupling /?, related to ry by = 1 — 2?7/7r. For the 
solution technique we wish to describe below, however, a rather less physical parameter will 
be especially important, for reasons which will be sketched shortly: the spectral parameter ^. 

To evaluate Z , we can first define the transfer matrix T, to be the s\im over the spins 
on a line of horizontal links of the product of their Boltzmann weights, for a given set of 




(1) 
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spins on the adjacent vertical links above and below. If we assign one multiindex a to the 

vertical links above the line in question, and another a' to the links below, then T can be 
thought of as a 2^ X 2^ matrix T^,, and, for periodic boundary conditions, the partition 
function Z is 

Z = Ttace[T^] . (2) 

Clearly, if we can diagonalise T, then we shall be able to compute Z. Even better, since the 
main interest is in the so-called thermodynamic limit when both N and M tend to infinity, 
we only need find the lowest-lying eigenvalues. Note that the mathematical structures found 
in this limit can also be obtained directly in the context of the continuum quantum field 
theory of the massless sine-Gordon model on a cylinder, using the constructions of [9-11]. 

Now comes the first moment where the concept of integrability arises: the local Boltz- 
mann weights for the six-vertex model are so defined that the transfer matrices taken at 
different values of the spectral parameter commute: 

[T(C),T(f')]=0. (3) 

(The deeper meaning of this property would take us too far afield, but from one point of 
view it is related to the existence of infinitely-many commuting conserved quantities.) If 
the matrices T(^) all commute, then they can be diagonalised simultaneously and we can 
study individual eigenvalues T(^) as functions of ^. There are various methods for finding 
these functions, but the most relevant here was discovered by Baxter in the early 1970s, in 
connection with his work on the more-complicated eight-vertex model. Since the transfer 
matrix is an entire function of ^ with ^-independent eigenvectors, the eigenvalues T(^) must 
also be entire. Baxter showed further that, for each eigenvalue T(^), there is an associated 
function Q{S,), also entire, such that the following TQ relation holds: 

nOQiO = e-2-PQ(g-2e) + e'-^PQiq^O (4) 

where the phase q depends on /3 as g = e'^'^ , and the extra parameter p allows for the pos- 
sibility of inserting a 'twist' into the periodic boundary conditions on the lattice. (Strictly 
speaking, for finite values of the lattice width N the TQ relation contains some extra fac- 
tors, but these go in the thermodynamic limit. We are following the conventions of [10], 
appropriate for this limit, save that there Q is denoted A for p ^ 0, and both it and T are 
taken to be entire functions of the squares of their arguments.) To see why this equation 
provides a useful constraint, we can argue as follows. Since Q is entire, it is, subject to 
suitable growth conditions, determined by the location of its zeroes. Suppose these zeroes 
are at {eo . . . e/c} ; then 

Q(o=Q(o)n(i-f) . (5) 

i=0 ^ 

Now T(^) is also entire, so the TQ relation (4) taken at ^ = ej implies that e~^'^''PQ{q~^ej) + 
e^^^PQ{q^ej) = 0. Rearranging and using the product form (5), 

/ 2 
TT I e» - g ej 

. „ V Ci — a~^e,- 




(6) 
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This gives us k+1 equations in k+1 unknowns, which determine the ej up to discrete ambi- 
guities corresponding to the fact that the transfer matrix has not one but many cigcnvahies. 
These equations coincide with the so-caUed Bethe ansatz equations (BAE) found in other, 
more direct, approaches to the problem. In this context, the Cj are sometimes called the 
Bethe roots. Note, once the Bcthc roots arc known, the eigenvalue T(^) is easily recovered 
by using (5) and then (4). Although the BAE (6) do not determine the Bethe roots uniquely, 
for the so-called ground state eigenvalue - the largest one, most important in the thermo- 
dynamic limit - this ambiguity is fixed by the fact that the Bethe roots for this eigenvalue 
are all located on the positive real axis. 

Over the last few decades, Bethe ansatz equations such as (6) have been studied exten- 
sively, and much is known about their solutions. For the moment the key thing to remember 
is the simple way that they follow from the TQ relation. If a similar functional equation 
can be found in another context, we can hope to establish a rather detailed connection with 
the already-existing body of knowledge on integrable models. 

3. Second prologue: PT^symmetric quantum mechanics 

We now make a change of tack, and turn our attention to some quantum-mechanical prob- 
lems first investigated by Bessis and Zinn- Justin in the early 1990s [3], and then successively 
generalised in [4,5,17,24,25]. For the current discussion it will suffice to consider the fol- 
lowing collection of non-Hermitian 'position space' Hamiltonians: 



where M is a positive real number. The interest is in the spectrum of Hm,i ■ those values of 
E such that the ordinary differential equation HM,ifp = Etp has a normalisable solution. By 

'normalisable', for Af < 2 we simply mean square-integrable on a contour running along the 
real axis, distorted just below the origin to avoid the singularity there when ^(/+1) ^ 0. For 
M > 2, the contour must be further distorted to ensure the correct analytic continuation, 
as explained in [4]. 

In figure 1 we show the spectrum as a function of M for two negative values of I. Note 
that M = 1 corresponds to the exactly-soluble simple harmonic oscillator with angular 
momentum, and that this point marks a profound change in the nature of the spectrum: for 
Af > 1 it is entirely real, while for M < 1 infinitely-many eigenvalues pair off and become 
complex. The reality of the spectrum for M = 3/2, / = was the subject of Bessis and 
Zinn- Justin's original conjecture [3], which was extended to all M > 1, I = by Bender and 
Boettcher in [4]. However a complete proof was surprisingly elusive, and has only recently 
been given, making essential use of the ODE/IM correspondence [24]. 

The nature of the transition to complex eigenvalues as M falls below 1 looks rather 
simple on figure la, but less so once figure lb has been examined. In fact, as shown in figure 
2, for Z > the connectivity pattern seen in figure la is completely reversed. Some further 
discussion of this behaviour can be found in [17,29], while the even-richer structure which 
emerges when a suitably-chosen inhomogeneous term is added to the potential is discussed 



As emphasised in [4,5], a key feature of the Hamiltonians (7) is that, while not Hermitian, 
they are unchanged by the combined action of parity P and time-reversal T - whence the 




(7) 



in [24,25]. 
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M M 
la: / = -0.025 lb: I = -0.001 



Figure 1. Eigenvalues of Ti-M^ii} = ^V") plotted as a function of M at i = —0.025 and I = —0.001. 
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2a: / = 2b: I = 0.001 



Figure 2. Eigenvalues of Ti,M,lip = Etp, plotted as a function of M at i = and I = 0.001. 

name 'PT-symmetric quantum mechanics. We hope to have demonstrated that the subject 
contains many interesting phenomena for the mathematical physicist to expfore; see [30-49] 
and references therein for some further devefopments in the area. However, from the point 
of view of this article, its main interest lies in its surprising links with the six-vertex model 
discussed in section 2 above, which we shall now describe. 
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4. Spectral determinants and the correspondence 

Rather than worrying about individual eigenvalues of (7), it turns out to be useful to com- 
bine them into a single function - a spectral determinant - and then exploit the analytic 
properties of this function. This approach was particularly advocated in the differential 
equation context by Sibuya [1] and Voros [2] . Here, influenced by parallels with the theory 
of integrable models, we shall use spectral determinants to couple the non-Hermitian prob- 
lem discussed so far with a related Hermitian problem, obtained from (7) by sending x ^ x/i 
and E ^ —E. The new problem is 

-^+-^" + ^)*(-)=^^*(-)> (8) 

and we choose as boundary conditions that the solution should vanish as x oo along the 
real axis, and behave 0. For ^Rel > —1/2, this problem is Hermitian. In the 

language of ordinary differential equations in the complex domain, it is sometimes called a 
'radial' problem, while (7) is called a 'lateral' problem. (Note, the radial problem can also 
be considered for 'Re I < —1/2, but is best then defined by analytic continuation in I.) 

Let {Ej} be the set of eigenvalues of (7), and let {cj} be the eigenvalues of (8). Then 
define two spectral determinants, as follows: 

T{E) = mf[(l + ^^ , (9) 

and 

Q{E)=Q{0)f[(l-^) . (10) 
i=o ^ 

Both products are convergent for M > 1, and define entire functions of E, and their zeroes 

(or, for T{E), their negatives) coincide with the eigenvalues of the corresponding spectral 
problem. For M < 1 convergence factors must be added, and it is more efficient to define 
the spectral determinants indirectly, via certain special 'Sibuya' solutions to (8), which are 
anyway needed to prove the key identity (11) below; see [1, 17] for more details. 

By considering the asymptotic behaviour of the Sibuya solutions in the complex plane it 
is not hard, following the arguments given in [17], to establish a Stokes relation, from which 
one can obtain the following functional equation 

T{E)Q{E) = a;-(2'+i)/2g(a;-2i;) + uj^-^^+^^^Qii^'^E) , (11) 

where u) = e"^/(-'^+^). This shows that the spectral problems (7) and (8) are related by 

much more than a simple change of variables and boundary conditions; and furthermore, 
that this relation is encoded precisely in a TQ relation of the sort which had previously 
arisen in the context of the six- vertex and sine- Gordon models. 



5. Consequences of the correspondence 

If we set = 1/(M+1) and p = (2Z-|-l)/(4M-|-4), then the two TQ relations (4) and 

(11) match perfectly. The remaining ambiguity in the solutions of the resulting Bethe 
ansatz equations is resolved once we recall that the problem (8) is Hermitian, and so all 
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of its eigenvalues are real. Given the remarks about the locations of Bethe roots at the 
end of section 2. this means that the spectral determinant Q{E) arising in the ordinary 
differential equation should be identified with the Q-function Q{i) for the ground state of 
the sine-Gordon model. (We must also match the normalisations of E and ^, since these 
are left undetermined by the TQ relations. This is easily achieved by comparing the leading 
asymptotic of the field theory Q-function given in [10] with that obtained by a simple WKB 
analysis of the ordinary differential equation - see [6, 17] for more details.) Once the Q- 
functions have been identified, so can be the T-functions, and we have our main punchline: 

The spectral determinant T for the lateral ('PT-symmetric) problem (7) is equal to 
the ground state, or vacuum, eigenvalue of the transfer matrix T of the massless 
sine-Gordon model on a cylinder. 

Together with the corresponding statement for Q, this provides a powerful tool to analyse 
the functions constructed in the context of integrable quantum field theory in [9-11], as 
the analytic properties of solutions of ordinary differential equations are relatively easy to 
control. (In this regard, we should mention that, for a complete proof of the correspondence, 
it currently seems best to argue via the so-called 'quantum Wronskian' relation, as in [6,14]. 
However the route we have adopted here is perhaps more intuitive, and it leads naturally on 
to the reality proof which we shall describe shortly.) There arc also less direct applications 
of the correspondence - for example, it yields a proof [14] of a duality relation in quantum 
Brownian motion, first proposed in [50]. In finishing this contribution, we mention two 
more, both on the ordinary differential equations side of the correspondence. 

The first concerns the reality properties of the PT-symmctric problems (7). We have 
already seen, in section 2, the argument which shows that substituting E = Cj into the TQ 
relation (11) leads to Bethe ansatz equations, which in the light of the correspondence we 
can reinterpret as coupling the different eigenvalues of the Hcrmitian problem (8). However, 
if we instead set E ^ —Ej, the same reasoning leads, in the region M > 1 for which the 
factorised form (10) apphes, to 



This equation couples the so-far mysterious eigenvalue Ej of the PT-symmetric problem (7) 
to the much better-controlled, and indeed real and positive, eigenvalues of the Hermitian 
problem (8). If we take the modulus^ of both sides, it is soon seen that the only way that the 
resulting equality can be achieved is for Ej to be real, as had previously been conjectured. 
Notice, the proof breaks down for M < 1, as has to be the case given the numerical findings 
of [4,17], illustrated in figures 1 and 2 above. More details, in a slightly more general setting 
than that described here, can be found in [24]; for a further generalisation of the method, 
see [51]. 

The second application was also found in [24], but is a little more specialised. We 

consider the radial problem at M = 3, and add an inhomogeneous term ax"^ to the potential. 
Although this was not done in the original paper, it will also be convenient to reparametrise 




(12) 
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the angular-momentum term by setting A = \/3(2Z-|-l) , so that the equation becomes 

Via the Bethe ansatz approach, it turns out that this problem has a relationship with a 
thirc^oidei ordinary differential equation: 

{D{g2-2)D{g,-l)D{go) + x^) <f> = ^ Ecj^ (14) 

where Dig) = {£ - f ), and 

50 = 1 + (a + V3A)/4 , gi = 1 + a/2 , g2 = 1 + {a - V3X)/4 . (15) 

This third-order equation is associated with SU{3) Bethe ansatz equations, as discussed 
in [18,21]. Furthermore, the third-order equation is symmetrical in {go, gi, g2}, a f^eature 
which is completely hidden in the original second-order equation. By playing with this 
symmetry, one can establish some novel spectral equivalences between different (second- 
order) radial problems, and also between these and certain lateral problems. More details 
and explicit formulae can be found in [24] , and here we simply add to that paper the remark 
that, when expressed in terms of the variables (a, A), the mappings turn out to act as certain 
2x2 matrices in the Wcyl group of SU{3). 

Other aspects of the ODE/IM correspondence are still being developed as this contri- 
bution is being written. It would be nice to have a working correspondence also for finite 
lattice systems, before the thermodynamic limit is taken, and some small progress in this 
direction will be reported in [52]. Another important question that imtil recently remained 
open concerns the other states in the quantum field theories, besides the ground state, and 
whether they can also be matched with differential equations. In a very recent paper [53], 
Bazhanov, Lukyanov and Zamolodchikov have answered this in the affirmative. The 'Q- 
potentials' that they construct are no longer real, even for the radial problems. This is not 
as surprising as it might seem - recall that it is only for the ground state that the Bethe 
roots, which correspond to the eigenvalues of the radial spectral problem, all lie on the real 
axis. For other states there are complex roots, and so even the radial problems can no longer 
be Hermitian. 

Finally, it would be valuable to have a more physical understanding of the relationship 
between integrable quantum field theories and ordinary differential equations. At the mo- 
ment this remains mysterious, but sufficiently-many examples of the phenomenon have now 
been collected that progress may not be too far off. 
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